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Introduction: Supernova explosions are the most violent and energetic events in the galaxy 
and have long been considered probable sources of Cosmic Rays (1). Recent shock accelera- 
tion models (2,3), treating the Cosmic Rays (CR's) as test particles in a prescribed 
Supernova Remnant (SNR) evolution, indeed indicate an approximate power law momentum 
distribution f SO urce^ p ^ p " a for the P art icles ultimately injected into the Interstellar 
Medium (ISM). This spectrum extends almost to the momentum p = 10 6 GeV/c, where the break 
in the observed spectrum occurs. The calculated power law index a ~ 4.2 agrees with that 
inferred for the galactic CR sources. The absolute CR intensity can however not be well 
determined in such a test particle approximation. 

It is important to know the intensity which results if the CR's are selfconsitently 
included in the dynamics of the SNR, and whether the predicted spectral characteristics 
are approximately conserved in such a nonlinear treatment. In the average ISM the CR 
energy density is comparable to the thermal and magnetic energy densities. Thus, in their 
sources the CR's must also be a significant dynamical element that manifests itself 
through its pressure as well as its high energy current (4,5). In fact the usual estimates 
for the observationally required CR production efficiency of SNR sources are too low by 
about a factor 4 which raises the overall production efficiency from the traditional 2 to 

4 percent of the average SNR energy input rate 
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1: The inferred source distribution 
f source foT ri 9 idit y dependent escape 
from the galaxy up to the cutoff imposed 
by the finite acceleration time in 
SNR's. 


but the source intensity f source , obtained 
from f obs by correction for rigidity dependent 
escape r e (p)^p" e must be replenished after a 
CR lifetime r g (p rz l GeV/c) of about 2-10 7 yrs. 
Taking £= 0,55 (6,7,8) and a cutoff momentum 
of 10^ GeV/c, the additional input rate corre- 
sponding to the hatched region in Fig. 1 is 
about a factor of 2.75 larger than the one 
corresponding to f Qbs , leading to a factor 
(4.75-4)/(4.75-4- 6) = 3.75 increase in the 
overall efficiency requirement. According to 
the test particle models, most of the accel- 
eration must take place prior to the phase in 
which the SNR radius is about a factor 4 below 
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its radius at particle release. Therefore the instantaneous early conversion efficiency 
must be about 4 times higher than the required average acceleration efficiency and corres- 
pond to 32 to 64 percent of the total energy input rate per supernova. Thus one must 
expect an extensive CR reaction on SNR dynamics. This is not only of interest for CR 
origin but also for the overall structure of the ISM in galaxies. 

This paper discusses simplified models for SNR evolution with CR's. Although ultimately 
this must be calculated numerically, such models are needed to clarify the numerical 
results. They are important for parameter studies and provide a first iteration and rough 
guidance for the interpretation of observations. Since the present theory only calculates 
macroscopic CR quantities we do not get spectral information about the particle distrib- 
ution, but information on the most important unknown, the overall efficiency of the 
acceleration process. 


2. Phases of SNR evolution: In the initial phase of a SNR the cooling ejecta sweep up 
external interstellar material until a reverse shock runs into the interior and heats it 
to high temperatures. Even though the reverse shock will certainly accelerate particles 
from the ejecta and may in fact push the maximum particle energy from a disk SNR easily to 
10^ GeV with interesting consequences for the chemical and isotopic composition in this 



Fig. 2: Density field in a 2-shell model 
for a SNR, where upstream ISM, shell, 
and interior parameters are denoted by 
the suffices 1,2, and 3 respectively. 

All mass is concentrated in the shell, 
all internal energy in the interior. 

Mass balance is then given by 


energy range, it is not included in the model 
presented here. We consider the subsequent 
phases where the outer shock is driven by a 
uniform hot interior generating a uniform 
dense shell with a thickness &R that is small 
compared to the remnant radius R(t) (Fig. 2). 
The local physical variables in the shell are 
assumed to have their postshock values. 

Such an approximation is by no means new. For 
the adiabatic phase it has been used e.g. by 
Chernyi (9) and it is standard for the so- 
called pressure modified snowplow phase after 
the onset of shell cooling. We apply it here 
for particle acceleration and backreaction on 
the gas dynamics in a hydrodynamic approxima- 
tion for the CR, or High Energy Component as 
the energy containing part of the CR spectrum 
below 10 6 GeV/c might be more appropriately 
called (10). 
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, where the total mass M g is assumed to be concentrated in the shell and dM^/dt ~T 3 
denotes the mass evaporation rate from cold clouds in the hot interior with temperature 
T y Momentum balance is that of the shell 
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where v 2 is the postshock mass velocity, and the total pressure p = p + p is the sum of 
gas pressure p g and CR pressure p Q . In the adiabatic phase where gas Pooling is unimpor- 
tant overall energy conservation - neglecting any other CR energy losses than adiabatic 
ones - reads as 
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Herey- 5/3, and 4/3 <y c < 5/3 are the effective adiabatic indices of thermal gas, and 
CR's, respectively. CR energy balance can be written as 
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, where F Q = v.p c .y c /(y c -l) - (4p c /4r)-i?/(y c -l) is the local CR energy flux density, H the 
mean CR diffusion coefficient (10), and we have assumed v(r,t) = v 2 (t).r/R(t) in rough 
correspondence to the selfsimilar Sedov solution. 

For the particle acceleration at the shock (Fig. 2) we assume a quasi-stationary state, 
i.e. steady fluxes of mass, momentum and total energy in the shock frame plus local 
Energetic Particle Hydrodynamics (11). Due to the finite acceleration time this is a far 
reaching assumption indeed. However, adiabatic cooling in the expanding SNR introduces a 
finite postshock gradient in CR pressure which couples the acceleration at the shock to 
the interior dynamics. Among other effects this is equivalent to an instantaneous cutoff 
of the accelerated spectrum at that particle momentum where the acceleration time becomes 
equal to the SNR lifetime. 

Radiative cooling of the shock heated gas sets in for T ^ 10 6 K. Unless cloud evaporation 
were so efficient that cooling set in first in the interior (12), cooling leads to a dense 
shell and an adiabatic interior: P g3 ~R" 5 and v 2 ~dR/dt. Of course the CR's do not cool 
radiatively. Due to the now lower shock speed some of the highest energy particles escape, 
but the rest still continues to be accelerated, driving the remnant. The higher energy 
particles should then be accelerated across the dense shell where they are scatterfree due 
to ion-neutral damping of hydromagnetic waves (13). Neglecting P g2 -P g3 across the thin 
shell, continuity of mass and momentum flux allows one to simply determine the pressure 
P C 2 of accelerated particles 


p c2 - • (dR/dt) 2 . p c1 ♦ <P gl - P g3 > 
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in terms of the shock velocity and overall gas pressure contrast. 
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3. Preliminary conclusions: (i) the CR’s take away internal energy from the gas with 
earlier onset of cooling, while on the other hand (ii) the CR’s cool much less due to 
their softer equation of state so that they can push the shell from the interior longer 
than the gas (iii) coupling with the interior introduces an effective time dependence of 
acceleration through a cutoff in the spectrum (iv) disregard of the sweep-up phase makes 
the model somewhat inconclusive regarding the relative contribution of gas and CR pressure 
during the earliest phases (v) After cooling the effective shock compression ratio in- 
creases, leading possibly to harder momentum spectra for older shocks with intriguing 
consequences for the radial (synchrotron) spectral index variation in SNR's. 
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